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ABSTRACT 

We fitted the optical to mid-infrared (MIR) spectral energy distributions (SEDs) of ^15000 type-I, 0.75 < 
z < 2, active galactic nuclei (AGNs) in an attempt to constrain the properties of the physical component re- 
sponsible for the rest-frame near-infrared (NIR) emission. We combine optical spectra from the Sloan Digital 
Sky Survey (SDSS) and MIR photometry from the preliminary data release of the Wide Infrared Survey Ex- 
plorer (WISE). The sample spans a large range of AGN properties: luminosity, black hole mass, and accretion 
rate. Our model has two components: a UV-optical continuum source and very hot, pure-graphite dust clouds. 
We present the luminosity of the hot-dust component and its covering factor, for all sources, and compare it 
with the intrinsic AGN properties. We find that the hot-dust component is essential to explain the (rest) NIR 
emission in almost all AGNs in our sample, and that it is consistent with clouds containing pure-graphite grains 
and located between the dust-free broad line region (BLR) and the "standard" torus. The covering factor of 
this component has a relatively narrow distribution around a peak value of ~0.13, and it correlates with the 
AGN bolometric luminosity. We suggest that there is no significant correlation with either black hole mass or 
normalized accretion rate. The fraction of hot-dust-poor AGNs in our sample is ~ 15 — 20%, consistent with 
previous studies. We do not find a dependence of this fraction on redshift or source luminosity. 

Subject headings: catalogs — galaxies: active — infrared: galaxies 



1. INTRODUCTION 

The unification scheme of active galactic nuclei (AGNs) re- 
quires an anisotropic obscuring structu re that surrounds the 
central accreting black hole (e.g., Krolik & Begermanl ll988t 
I AntonuccH 1 1 9931) . In this picture, the bulk of the radiation 
from the central engine is absorbed by the obscuring structure 
(commonly referred to as the torus) and re-emitted mainly in 
mid-infrared (MIR) wavelengths. 

A component of very hot dust at the inn ermost edge of the 
torus has been suggest ed in the past (e.g., iNeugebau er et al.l 
119871: Barvainisl 11987b . and in recent years is increasingly 
being supported by observations. Reverberation measure- 
ments of nearby AGNs suggest that the near infrared (NIR) 
emission in these sources is dominated by thermal radiation 
from hot dust very close to the central source (few tens o f 
light days; e.g., Minezaki "et alJl2004tlSuganuma et"aT|[ 2006). 
Other studies fitted the NIR-MIR spectral energy distributions 
(SEDs) of AGNs using a blackbody spectrum to represent 
emission from hot dus t in the inner region of the toru s (e.g., 
iKishimoto etal] [20071: iRiffel et alj [20091 : iMoretalJ 120091 
hereaf ter M09; (Deo et al.ll201 ll) .~~More recently lLandt et all 
( 1201 11) found similar results by fitting only the optical-NIR 
SED of 23 AGNs. The modeled temperature of this compo- 
nent is found to be high, £ 1200K, regardless of the AGN 
luminosity and consistent with pure-graphite dust composi- 
tion (M09). Several studies have shown that the luminosity of 
the NIR e xcess emission correl ates with that of the central en- 
gine (e.g., Gallagher et al. 2007, and references therein) with 
a slope close to unity. However, it is yet unclear whether this 
hot-dust component is related to other AGN properties such 
as its mass (Mbh) or normalized accretion rate (L/Lndd)- 

Although hot dust seems to be a common feature of AGNs, 
several recent studies have suggested that a certain fraction 
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of the AGN population lacks such a component. Jian g et alj 
(2010) found two z~6 QSOs without any detectable emission 
from hot dust (using Spitzer MI R photometry). T hese were 
dubbed "hot-dust free" Q SOs. lHao et al.1 (120101) and more 
recently lHao et all (1201 ll) found a sizable amount of type- 
I AGNs with unusually weak NIR emission in several large 
samples of type-I AGNs. These were dubbed "hot-dust-poor" 
(HDP; hereafter we adopt this notation) AGNs. These au- 
thors found that the fraction of HDP AGNs increases with 
redshift from 6% at low redshift (z < 2) to 20% at higher red- 
shift (2 < z < 3.5). We note that all the studies mentioned 
here are limited to relatively small samples, and use highly 
simplified emission models. 

In this letter we use the recently published preliminary 
data release of the Wide Infrared Survey Explorer (WISE; 
Wright et al. 2010), together with the sevent h data release of 
the S loan Digital Sky Survey (SDSS/DR7; lAbazaiian etal] 
[20091) to construct - 15000 UV to NIR SEDs of type-I AGNs. 
We apply spectral decomposition using novel models to de- 
duce the properties of the hot dust, and test these properties 
against different intrinsic AGN properties. 

2. SAMPLE SELECTION AND SPECTRAL DECOMPOSITION 

We aim to study a large and uniformly selected sample of 
type-I AGNs, for which Lbob -^bh an d ^/^Edd can be reliably 
measured, and the WISE bands (at ~ 3.4,4.6, & 11.6 jUm) 
cover the rest-frame wavelength range of 1 to 5 /xm, where the 
emission originating from the hot dust is expected to peak. At 
longer wavelengths (i.e., those covered by the ~ 22 jim WISE 
band), the emission from the "standard" torus becomes sig- 
nificant (see Fig. 4 in M09). Thus, we select from the SDSS 
only high confidence "QSO" sources (i.e. zconf > 0.7) at 
0.75 < z < 2. This query returned 20231 sources within the 
area covered in the WISE preliminary DR. Cross-matching 
this sample with the WISE catalog, using a 3" search radius, 
yielded 19116 matched sources. The remaining 1115 sources 
(~ 5.5%) probably lie below the flux limit defined for the 
WISE preliminary DR, of ~ 0.08 mJy in the 3.4 jum band. 
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These sources have a similar range of Lb i to that of the whole 
parent sample, and are further discussed in $3] Thus, our anal- 
ysis of the WISE data provides useful information regarding 
~ 94.5% of optically selected, 0.75 < z < 2, type-I AGN (see 
more criteria below). 

We further filtered the resultant cross-matched catalog to 
include only the 17920 sources which have S/N > 3 in each 
of the ~3.4, 4.6 & 11.6 jUm WISE bands. Our analysis ig- 
nores the ~ 22 jim band. The WISE cataloged magnitudes 
and associated uncertainties were translated to f v through the 
published ("iso") WISE zero points. Q The SDSS spectra 
of all these QSOs were modeled by a comprehensive pro- 
cedure that is aimed at fitting the broad Mg II A2798 emis- 
sion line and the adjacent emissio n complexes. The proce - 
dure is identical to that used in Trakhtenb rot et al.l (120111). 
and s imilar to other studies (e.g.. lShen et al.ll2.Q08l; iFine et alj 
2008). Most importantly, it allows a robust determination 
of FWHM(Mg II) and of the monochromatic luminosity at 
(rest-frame) 3000A (A L A [3000 A], hereafter L 30u o)- Further 
1571 sources, for which the Mgll fitting procedure failed 
to provide high-confidence results, were excluded. We also 
excluded 421 sources with log (Z^ooo/ergs -1 ) < 45, which 
might contain some host-galaxy contribution to their rest- 
frame NIR SED, and often have SDSS spectra with low 
S/N. Finally, our sample consists of 15928 QSOs with 45 < 
log (Z^ooo/ergs -1 ) < 47.25, which correspond to about 78% 
of all SDSS selected 0.75 < z < 2, type-I AGN covered by 
WISE observati ons. From these measur ements, we deduce 
Mbh using the McLure & Dunlopl (120041) prescription. Lbol 
is calculated from L3000 by applying a luminosity-dependant 
bolometric correction, which will be discussed in full de- 
tail in a forthcoming publication. Here we only note that 
these correction factors were independently calibrated using 
a large sample of 0.5 < z < 0.75 QSOs, f or wh ich Lbol could 
be reliably derived by the iMarconi et al.l (12004b prescription. 
The bolometric corrections applied to the SDSSAVISE sam- 
ple range between about 3, for the most luminous sources, 
and 4.2, for the faintest. These v alues are lower, by a f ac- 
tor of ~ 1 . 5 - 1. 7, than those of iRichards et all (120061) or 
lElvis et alJ 0994), since they avoid the double-counting of 
the re-processed IR emission (see the detailed discussion in 
IMarconi et al.l20 04) . The resultant range in log (Lbol / erg s~ 1 ) 
is 45.6 —47.7. We estimated the accretion rate as L/Lem = 

Lbol/ (1.5 xl0 38 M BH ). 

The SDSSAVISE SED of each source in our sample is 
fitted to constrain the properties of the hot-dust compo- 
nent that dominates the NIR wavelength range. To focus 
on continuum emission, every SDSS spectrum was sampled 
at several continuum bands, around the (rest frame) wave- 
lengths: 1450 - 1475, 2150 - 2200, 3030 - 3100, 4150 - 
4250, 5080 - 5120, 5600 - 5750, &6100 - 6250 A (e.g., 
IVanden Berk et al] 12001). Our models include two differ- 
ent components: a source of UV-optical continuum emission 
(e.g., an accretion disk) and very hot-dust clouds with pure 
graphite composition. 

To model the central-source continuum emission we as- 
sume a power-law function, L v V a . This component has 
two free parameters, the power-law index and the normaliza- 
tion factor. The power-law index, a, is changed in steps of 
0.1 between -1.5 to 1. The second component represents a 
collection of dusty clouds of gas with a pure-graphite grain 

2 See: http://wise2.ipac.caltech.edu/docs/release/prelim/expsup/figures/sec4_3gi 
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FIG. 1.— Best fit to the SED of SDSS J134154.22+005948.6 (z = 1.715) 
using a combination of power-law (red line) and hot, pure-graphite dust (blue 
line) components. The original SDSS spectrum is shown in gray and the 
continuum bands selected for the fitting procedure are highlighted in black. 
The best fit model is shown in red. The quality of the fit is demonstrated by 
the ratio between the data and the model (bottom). 

compositio n. The SED o f this c omponent is taken from the 
models of Mo r & Netzerl (1201 ll) . which provide a full dis- 
cussion of all model assumptions and parameters. In short, 
this study adopts the "cloud model" of the broad-line region 
(BLR) and considers the graphite-containing gas to be the ex- 
tension of the BLR. The model assumes gas composition of 2 
Z Q , standard ISM-type dust depletion, gas density in the range 
log(n/cm~ 3 ) = 9.3 —9.8 and column density in the range 
log = 22.3 — 22.7 at the graphite sublimation ra- 

dius. While the local emission of the graphite dust depends 
only on the grain properties, the dust temperature inside the 
cloud varies by a large factor because of the local grain opac- 
ity. Thus a single cloud spectrum appears as a combination 
of many modified blackbodies. The results of the calculations 
provide a grid of SEDs that are used in the fitting procedure. 
For this component there are two free fitting parameters, the 
distance to the cloud (given a source luminosity) and a nor- 
malization factor which determines its luminosity (hereafter 
Lhd)- For example, the range of distances for an AGN with 
log (Lbd/ergs -1 ) = 46.9 is 1.3 to 3.5 pc. The fitting proce- 
dure uses a standard % 2 minimization to determine the best 
fit combination of power-law and hot-dust model. Figure Q] 
demonstrates the fitting procedure for a representative case. 
The hot pure-graphite dust component dominates the SED be- 
tween ^2 and 5 fim. The main caveat of the fitting procedure 
is the (unknown) contribu tion to the SEP fro m cooler dust in 
the torus. As shown in lMor & Netzerl (1201 ll) . neglecting this 
component may lead to an overestimation of Lhd by merely 
~ 10 — 20%. In the following analysis, we do not correct for 
this systematic uncertainty. 

The chosen redshift range, and thus typical AGN luminosi- 
ties, securely omits any host galaxy dominated SEDs from our 
sample. To minimize the effect of the (unknown) host-related 
extinction, we further limit our sample to include only sources 
for which the best-fit slope of the power-law component sat- 
isfies a > —0.8. This choice is comparable to those made in 
t4 (S(fier studies of the optical to NIR SED of AGNs. For exam- 
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FIG. 2. — Hot-dust luminosity vs. AGN bolometric luminosity. Results 
from the complete WISE sample studied in this work are shown as gray sym- 
bols. Black symbols highlight the sources with a > —0.4. The QUEST 
sample of PG QSOs are shown in red. The dashed lines represent Lhd x A>ol 
(i.e. constant CFhdX with normalization of 0.05 and 0.35. 



pie, the lVanden Berk et all (120011) composite has a ^ -0.44 
over the relevant wavelength range, while lHao et al.l (l20lob 
show that a rang es between about -0.8 and 0.6, based on the 
Elvis et all (119941) SEDs. In any case, omitting this criterion 
adds only a small fraction of sources and does not signifi- 
cantly change our main results. A more comprehensive anal- 
ysis of the extinction is beyond the scope of this Letter. The 
removal of these sources leaves 15077 sources that will be 
used in the following analysis. 

3. RESULTS AND DISCUSSION 

The most important parameter of the hot-dust component, 
determined by the fitting procedure, is Lhd- Figure |2] shows 
a prominent correlation between Lhd and Lb i- We also high- 
light a sub-sample of sources with extremely blue continua, 
which are defined as those with a > —0.4, and correspond 
to the bluest third of the sample (following Gallagh er et al.l 
2007). These sources are presumably less affected by extinc- 
tion. The distribution of these blue sources in the Lhd — Lb i 
plane is similar to that of the entire sample, suggesting that 
host extinction plays only a minor role in our analysis. Re- 
sults for the QUEST sample (PG QSOs studied in M09; red 
squares in Fig. [2j were calculated by employing the M09 pro- 
cedure (using the pure-graphite models) to the ^2 to 35 jJ.m 
SEDs. 

Assuming that the rest-frame NIR SED is due to repro- 
cessed AGN radiation by the hot-dust clouds, we can deduce 
the covering factor (CFhd) of the ce ntral source by these 
clouds, defined by CF H d = L H D/Lb i (Mai olino et al.ll2007T 
M09). Fig. |2] clearly indicates that CFhd of the entire popula- 
tion spans a relatively limited range. Indeed, the distribution 
of CFhd (Figure|3]l appears to be fairly concentrated around a 
peak value of 0.13 (median value is 0.1 with a standard devi- 
ation of 0.4 dex). This value is smaller than the one found for 
the lower luminosity, lower redshift, QUEST sample (^0.23). 
The distribution of CFhd is asymmetric with a clear excess of 
sources with low CFhd- These sources may represent a differ- 
ent subset of the population with a very small or no hot-dust 
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FIG . 3 . — Distribution of the hot-dust covering factor. Values of the current 
study (blue) are generally smaller than those found for the QUEST sample 
(red) and may be the result of the different luminosity range of the two sam- 
ples (see Figure[4}. The CFhd distribution is asymmetric with a clear excess 
of low CFhd sources. The peak of the distribution is at ~0.13 and the median 
value is 0. 1 . Sources with CFhd below a certain luminosity-dependant value 
may represent the so-called hot-dust-poor or dust-free AGNs. 



component. We note that for a small fraction of the sample 
(~ 1 .7%, 258 sources) the fit does not require a hot-dust com- 
ponent thus resulting in Lhd = and CFhd = 0. We add these 
sources to the sub-group of HDP AGNs discussed below. 

Fig- E] presents a clear anti-correlation between CFhd an d 
Lboi- This trend is further confirmed by both Pearson's 
and Spearman's rank correlation tests (p value <C 0.01 for 
both tests). In particular, the typical CFhd decreases from 
~ 0.2 for sources with log (Lb i/ergs~') ~ 45.6 to ~ 0.09 

for sources with log (Lboi/ergs -1 ) ~ 46.6. We verified that 
a similar relation is recovered even if one calculates Lboi 
with the commonly used, uniform (and overesti mated; see 
bolometric correction o f 5.15 or 5.62 (e.g., lElvis et ail 
119941: iRichards et al.l 12006 ). Several earlier studies sug- 
geste d a similar trend (e.g ., IWang et al.ll2005l iMai olino et al. 
2007; Treister et al. 2008), however these were based on the 



total MIR emission which translates to the covering fac- 
tor of the entir e dusty structure, and not just the hot dust. 
Gallagher et al.l d2007l) suggest that this CF H d — Lb i anti- 
correlation may be a manifestation of dust extinction. We 
therefore verified that this trend persists for the sub-sample of 
blue sources, defined above. We further checked whether the 
CFHD-Lbol relation is due to the higher abundance of more lu- 
minous AGNs towards higher redshifts, using a sub-sample of 
3837 sources which reside in a narrow luminosity range (0.2 
dex) around log (Lboi/ergs -1 ) ~ 46.3. Although this sub- 
sample spans almost the entire redshift range of the original 
sample, no correlation is found between CFhd and redshift. 

The physical mechanism responsible for the decrease of 
covering factor with Lb i is still un determined. On e possi- 
bility is a "receding torus" scenario (lLawrencelll991l) . where 
higher luminosity implies larger dust sublimation distance, 
and hence an obscuring structure that is located farther away 
from the center. In this scenario, however, the geometry of the 
hot-dust clouds must be toroidal and have a constant height. 
Another possibility is that the CFHD-Lbol anti-correlation is 
analogous to the anti-correlation found between the equiva- 
lent width of different BLR lines (H/3 and C IV) and AGN 
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FIG. 4. — Hot-dust covering factor vs. AGN properties. Symbols as in Fig. [2] Left: CFhd decreases with increasing Lb i- Magenta markers represent the mean 
SED of Elvis et aH 11994) . scaled to log (Lboi/ergs ') = 46 (cross); the three examples of hot-dust-poor objects of Hao et al. (2010, triangles) and the two z — 6 
dust-free QSOs of Jiang et al. (2010, plus signs). The correlations with Mbh and ^/^Edd (middle and right panels) are much less significant. In all panels, dashed 
black lines represent the peak of the CFhd distribution and the 3o~-equivalent range - see text 



luminosity (e.g. jNetzer et al. 2004; B askin & Laor1l2005l) . 

The correlations between CFhd and Mbh and L/Ledd 
(Fig. |U middle and right panels) are much less significant. 
We interpret these as the result of the marginal dependence of 
both Mbh and ^/^Edd on the source luminosity. This is further 
confirmed by the fact that these marginal correlations disap- 
pear completely in the log (Lbd/ergs ) — 46.3 sub-sample 
mentioned above. However, the large scatter makes it difficult 
to securely determine which correlation is more fundamental. 
We suggest that the apparent CFhd is sensitive to the imme- 
diate irradiating flux from the central source, but does not de- 
pend on the evolutionary stage of the accreting super-massive 
black hole (SMBH). 

For a certain fraction of the sources, the hot-dust a compo- 
nent is much less prominent. First, as noted above, ~ 1 .7% of 
the sources have Lhd = and thus CFhd = 0. These sources 
spread the entire range of Lbol- Second, there is a non negligi- 
ble population of sources with significantly lower CFhd val- 
ues than the main locus of sources (Fig. There is no corre- 
sponding population of high-CFnD outliers. Fig. |4] illustrates 
this using the boundaries of the 99th percentiles of the CFhd 
distribution. The percentiles are calculated by assuming that 
the CFhd distribution should be symmetrical around the peak 
and mirroring the high-CFnD side of the distribution. This is 
done in each (0.2 dex) luminosity bin separately. We suggest 
that all the points which lie below the lower dashed line in 
Fig. |4]can be regarded as HDP AGNs. Thus, we are able to 
provide an HDP criterion which depends on Lbol- For exam- 
ple, for log (Lboi/ergs -1 ) = 46 all sources with CFhd ^ 0.06 
may be considered as HDP AGNs. 

The total fraction of HDP AGNs in our sample is ~ 16.2% 
(2439 sources out of 15077, including the 258 sources with 
CFhd = 0). This number i s in rough agreement with that re- 
ported by lHao et all (120101) and lHaoet aT] (120111) . We do not 
find a significant dependence of this fraction with either lumi- 
nosity or redshift. We stress that using a HDP criterion that 
is independent of AGN luminosity can result in false corre- 
lations of the fraction of HDP AGNs with either luminosity 
or redshift. The latter may be a consequence of selection bi- 
ases of high luminosity sources towards higher redshifts. As 



mentioned in $2] 1115 sources (~ 5.5% of the total SDSS 
detections) probably lie below the flux limit defined for the 
preliminary DR of WISE. These sources may also represent 
HDP AGNs, thus should be added to the total fraction men- 
tioned above (~ 16.2%). 

As mentioned in Sjl] Jia ng et all (12010) suggested that a sig- 
nificant fraction of z — 6 AGNs are hot-dust free, i.e. lacking 
any emission from a hot, dusty component. Moreover, these 
authors claim that such systems are not observed at lower red- 
shifts. Taken at face value, our CFhd = sources can be re- 
garded as the z < 2 counterparts of the z — 6 hot-dust free 
AGNs. These are ind eed rare at low red shifts. However, the 
hot-dust free AGNs of lJiang et al.l d2010) have only upper lim- 
its in the (rest-frame) IR bands, which place them in our HDP 
AGN regime (see Fig.|4]i. Thus, our sample might offer many 
m ore z < 2 sources th at correspond to the hot-dust free notion 
of lJiang et a l. (2010), from both the WISE-undetected sub- 
sample (~ 5.5%) and the HDP AGN sub-sample (~ 16.2%). 
We also note that variability between the SDSS and WISE 
observational epochs would not effect the estimation of the 
fraction of HDP AGNs for large samples. It may, however, 
account for some of the scatter observed in the symmetric part 
of the CFhd distribution. 

We have demonstrated that the rest-frame NIR emission of 
the vast majority (~ 80%) of type-I AGNs in our unprece- 
dentedly large sample can be explained by emission from hot, 
pure-graphite dust clouds. These clouds reside at distances 
greater by a factor of about 3 — 10 than the BLR clouds, and 
have a typical CFhd of ~ 0.13. The lack of a significant corre- 
la tion between CFhd and Mbh stands in contrast to the result 
of lJiang et alj (|2010), who suggested that such a relation only 
exists at very high redshifts (z — 6), where dusty structures 
around AGNs were not yet fully developed. Interestingly, we 
do not find any significant correlation between CFhd and ei- 
ther L/Lgdd or redshift, despite the fact that our sample effec- 
tively covers a wide range in these properties (0.05 < L/Ledd 
and 0.75 < z < 2). We suggest that previous reports of a 
decrease in CFhd, and an increase in the fraction of HDP 
AGNs, with increasing redshift, may be due to a combination 
of the selection of high-luminosity sources at high redshifts 
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and the CFnD-Lboi anti-correlation reported here. However, 
the small overlap in Lbol between the low-z QUEST sample 
and the SDSSAVISE sample, and the limited range of Lbol f° r 
SDSS/WISE sources, at a chosen redshift, inhibits us from 
testing this issue further. 

To fully understand the evolution of hot dusty structures 
around accreting SMBHs, all these trends should be tested in 
large samples of high redshift type-I quasars, fo r which Lbol, 
Mbh a nd L/Lfm were rel i ably measured (e.g. , iNetzer et al.l 
12001 IWillott et al.1 l20l(it ITrakhtenbrot etaLl UmW . As 
demonstrated in this Letter, WISE data provide the optimal 
way to preform such studies, owing to its wide areal cover- 
age and relatively deep flux limit. In addition, WISE covers 
the wavelength range where hot-dust emission dominates the 
SED of AGNs, thus providing more direct probes of this emis- 



sion at high-redshift sources than NIR data. 
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